Variation and plasticity in life-history traits and fitness of wild Arabidopsis thaliana populations are not related to their genotypic and ecological diversity.

BACKGROUND: Despite its implications for population dynamics and evolution, the relationship between genetic and phenotypic variation in wild populations remains unclear. Here, we estimated variation and plasticity in life-history traits and fitness of the annual plant Arabidopsis thaliana in two common garden experiments that differed in environmental conditions. We used up to 306 maternal inbred lines from six Iberian populations characterized by low and high genotypic (based on whole-genome sequences) and ecological (vegetation type) diversity. RESULTS: Low and high genotypic and ecological diversity was found in edge and core Iberian environments, respectively. Given that selection is expected to be stronger in edge environments and that ecological diversity may enhance both phenotypic variation and plasticity, we expected genotypic diversity to be positively associated with phenotypic variation and plasticity. However, maternal lines, irrespective of the genotypic and ecological diversity of their population of origin, exhibited a substantial amount of phenotypic variation and plasticity for all traits. Furthermore, all populations harbored maternal lines with canalization (robustness) or sensitivity in response to harsher environmental conditions in one of the two experiments. CONCLUSIONS: Overall, we conclude that the environmental attributes of each population probably determine their genotypic diversity, but all populations maintain substantial phenotypic variation and plasticity for all traits, which represents an asset to endure in changing environments.

Multilocus pathogenic variants contribute to intrafamilial clinical heterogeneity: a retrospective study of sibling pairs with neurodevelopmental disorders.

BACKGROUND: Multilocus pathogenic variants (MPVs) are genetic changes that affect multiple gene loci or regions of the genome, collectively leading to multiple molecular diagnoses. MPVs may also contribute to intrafamilial phenotypic variability between affected individuals within a nuclear family. In this study, we aim to gain further insights into the influence of MPVs on a disease manifestation in individual research subjects and explore the complexities of the human genome within a familial context. METHODS: We conducted a systematic reanalysis of exome sequencing data and runs of homozygosity (ROH) regions of 47 sibling pairs previously diagnosed with various neurodevelopmental disorders (NDD). RESULTS: We found siblings with MPVs driven by long ROH regions in 8.5% of families (4/47). The patients with MPVs exhibited significantly higher FROH values (p-value = 1.4e-2) and larger total ROH length (p-value = 1.8e-2). Long ROH regions mainly contribute to this pattern; the siblings with MPVs have a larger total size of long ROH regions than their siblings in all families (p-value = 6.9e-3). Whereas the short ROH regions in the siblings with MPVs are lower in total size compared to their sibling pairs with single locus pathogenic variants (p-value = 0.029), and there are no statistically significant differences in medium ROH regions between sibling pairs (p-value = 0.52). CONCLUSION: This study sheds light on the significance of considering MPVs in families with affected sibling pairs and the role of ROH as an adjuvant tool in explaining clinical variability within families. Identifying individuals carrying MPVs may have implications for disease management, identification of possible disease risks to different family members, genetic counseling and exploring personalized treatment approaches.

How do stochastic processes and genetic threshold effects explain incomplete penetrance and inform causal disease mechanisms?

Incomplete penetrance is the rule rather than the exception in Mendelian disease. In syndromic monogenic disorders, phenotypic variability can be viewed as the combination of incomplete penetrance for each of multiple independent clinical features. Within genetically identical individuals, such as isogenic model organisms, stochastic variation at molecular and cellular levels is the primary cause of incomplete penetrance according to a genetic threshold model. By defining specific probability distributions of causal biological readouts and genetic liability values, stochasticity and incomplete penetrance provide information about threshold values in biological systems. Ascertainment of threshold values has been achieved by simultaneous scoring of relatively simple phenotypes and quantitation of molecular readouts at the level of single cells. However, this is much more challenging for complex morphological phenotypes using experimental and reductionist approaches alone, where cause and effect are separated temporally and across multiple biological modes and scales. Here I consider how causal inference, which integrates observational data with high confidence causal models, might be used to quantify the relative contribution of different sources of stochastic variation to phenotypic diversity. Collectively, these approaches could inform disease mechanisms, improve predictions of clinical outcomes and prioritize gene therapy targets across modes and scales of gene function. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Predator selection on phenotypic variability of cryptic and aposematic moths.

Natural selection generally favours phenotypic variability in camouflaged organisms, whereas aposematic organisms are expected to evolve a more uniform warning coloration. However, no comprehensive analysis of the phenotypic consequences of predator selection in aposematic and cryptic species exists. Using state-of-the-art image analysis, we examine 2800 wing images of 82 moth species accessed via three online museum databases. We test whether anti-predator strategy (i.e., camouflage or aposematism) explains intraspecific variation in wing colour and pattern across northern hemisphere moths. In addition, we test two mutually non-exclusive, ecological hypotheses to explain variation in colour pattern: diel-activity or dietary-niche. In this work, taking into account phylogenetic relationships, moth phenotypic variability is best explained by anti-predator strategy with camouflaged moths being more variable in wing patterning than aposematic species.

Temporal phenotypic variation of spinach root traits and its relation to shoot performance.

The root system is important for the growth and development of spinach. To reveal the temporal variability of the spinach root system, root traits of 40 spinach accessions were measured at three imaging times (20, 30, and 43 days after transplanting) in this study using a non-destructive and non-invasive root analysis system. Results showed that five root traits were reliably measured by this system (RootViz FS), and two of which were highly correlated with manually measured traits. Root traits had higher variations than shoot traits among spinach accessions, and the trait of mean growth rate of total root length had the largest coefficients of variation across the three imaging times. During the early stage, only tap root length was weakly correlated with shoot traits (plant height, leaf width, and object area (equivalent to plant surface area)), whereas in the third imaging, root fresh weight, total root length, and root area were strongly correlated with shoot biomass-related traits. Five root traits (total root length, tap root length, total root area, root tissue density, and maximal root width) showed high variations with coefficients of variation values (CV ≥ 0.3, except maximal root width) and high heritability (H2 > 0.6) among the three stages. The 40 spinach accessions were classified into five subgroups with different growth dynamics of the primary and lateral roots by cluster analysis. Our results demonstrated the potential of in-situ phenotyping to assess dynamic root growth in spinach and provide new perspectives for biomass breeding based on root system ideotypes.

Phenotypic variation in biting behavior associated with differences in expression of olfactory genes in the vector mosquito, Aedes albopictus (Diptera: Culicidae).

We evaluated miRNA and mRNA expression differences in head tissues between avid-biting vs. reluctant-biting Aedes albopictus (Skuse) females from a single population over a 20-min timescale. We found no differences in miRNA expression between avid vs. reluctant biters, indicating that translational modulation of blood-feeding behavior occurs on a longer timescale than mRNA transcription. In contrast, we detected 19 differentially expressed mRNAs. Of the 19 differentially expressed genes at the mRNA level between avid-biting vs. reluctant-biting A. albopictus, 9 are implicated in olfaction, consistent with the well-documented role of olfaction in mosquito host-seeking. Additionally, several of the genes that we identified as differentially expressed in association with phenotypic variation in biting behavior share similar functions with or are inferred orthologues of, genes associated with evolutionary variation in biting behaviors of Wyeomyia smithii (Coq.) and Culex pipiens (Lin.). A future goal is to determine whether these genes are involved in the evolutionary transition from a biting to a non-biting life history.

Phylogeography and phenotypic wing shape variation in a damselfly across populations in Europe.

BACKGROUND: Describing geographical variation in morphology of organisms in combination with data on genetic differentiation and biogeography can provide important information on how natural selection shapes such variation. Here we study genetic structure using ddRAD seq and wing shape variation using geometric morphometrics in 14 populations of the damselfly Lestes sponsa along its latitudinal range in Europe. RESULTS: The genetic analysis showed a significant, yet relatively weak population structure with high genetic heterozygosity and low inbreeding coefficients, indicating that neutral processes contributed very little to the observed wing shape differences. The genetic analysis also showed that some regions of the genome (about 10%) are putatively shaped by selection. The phylogenetic analysis showed that the Spanish and French populations were the ancestral ones with northern Swedish and Finnish populations being the most derived ones. We found that wing shape differed significantly among populations and showed a significant quadratic (but weak) relationship with latitude. This latitudinal relationship was largely attributed to allometric effects of wing size, but non-allometric variation also explained a portion of this relationship. However, wing shape showed no phylogenetic signal suggesting that lineage-specific variation did not contribute to the variation along the latitudinal gradient. In contrast, wing size, which is correlated with body size in L. sponsa, had a strong negative correlation with latitude. CONCLUSION: Our results suggest a relatively weak population structure among the sampled populations across Europe, but a clear differentiation between south and north populations. The observed geographic phenotypic variation in wing shape may have been affected by different local selection pressures or environmental effects.

Non-genetic phenotypic variability affects populations and communities in protist microcosms.

Intraspecific trait variation (ITV), potentially driven by genetic and non-genetic mechanisms, can underlie variability in resource acquisition, individual fitness and ecological interactions. Impacts of ITV at higher levels of biological organizations are hence likely, but up-scaling our knowledge about ITV importance to communities and comparing its relative effects at population and community levels has rarely been investigated. Here, we tested the effects of genetic and non-genetic ITV on morphological traits in microcosms of protist communities by contrasting the effects of strains showing different ITV levels (i.e. trait averages and variance) on population growth, community composition and biomass production. We found that genetic and non-genetic ITV can lead to different effects on populations and communities across several generations. Furthermore, the effects of ITV declined across levels of biological organization: ITV directly altered population performance, with cascading but indirect consequences for community composition and biomass productivity. Overall, these results show that the drivers of ITV can have distinct effects on populations and communities, with cascading impacts on higher levels of biological organization that might mediate biodiversity-ecosystem functioning relationships.

Characterization of Phenotypic Variability in Becker Muscular Dystrophy for Clinical Practice and Towards Trial Readiness: A Two-Years Follow up Study.

Background: Becker muscular dystrophy (BMD) is a dystrophinopathy due to in-frame mutations in the dystrophin gene (DMD) which determines a reduction of dystrophin at muscle level. BMD has a wide spectrum of clinical variability with different degrees of disability. Studies of natural history are needed also in view of up-coming clinical trials. Objectives: From an initial cohort of 32 BMD adult subjects, we present a detailed phenotypic characterization of 28 patients, then providing a description of their clinical natural history over the course of 12 months for 18 and 24 months for 13 of them. Methods: Each patient has been genetically characterized. Baseline, and 1-year and 2 years assessments included North Star Ambulatory Assessment (NSAA), timed function tests (time to climb and descend four stairs), 6-minute walk test (6MWT), Walton and Gardner-Medwin Scale and Medical Research Council (MRC) scale. Muscle magnetic resonance imaging (MRI) was acquired at baseline and in a subgroup of 9 patients after 24 months. Data on cardiac function (electrocardiogram, echocardiogram, and cardiac MRI) were also collected. Results and conclusions: Among the clinical heterogeneity, a more severe involvement is often observed in patients with 45-X del, with a disease progression over two years. The 6MWT appears sensitive to detect modification from baseline during follow up while no variation was observed by MRC testing. Muscle MRI of the lower limbs correlates with clinical parameters.Our study further highlights how the phenotypic variability of BMD adult patients makes it difficult to describe an uniform course and substantiates the need to identify predictive parameters and biomarkers to stratify patients.

Detecting genetic effects on phenotype variability to capture gene-by-environment interactions: a systematic method comparison.

Genetically associated phenotypic variability has been widely observed across organisms and traits, including in humans. Both gene-gene and gene-environment interactions can lead to an increase in genetically associated phenotypic variability. Therefore, detecting the underlying genetic variants, or variance Quantitative Trait Loci (vQTLs), can provide novel insights into complex traits. Established approaches to detect vQTLs apply different methodologies from variance-only approaches to mean-variance joint tests, but a comprehensive comparison of these methods is lacking. Here, we review available methods to detect vQTLs in humans, carry out a simulation study to assess their performance under different biological scenarios of gene-environment interactions, and apply the optimal approaches for vQTL identification to gene expression data. Overall, with a minor allele frequency (MAF) of less than 0.2, the squared residual value linear model (SVLM) and the deviation regression model (DRM) are optimal when the data follow normal and non-normal distributions, respectively. In addition, the Brown-Forsythe (BF) test is one of the optimal methods when the MAF is 0.2 or larger, irrespective of phenotype distribution. Additionally, a larger sample size and more balanced sample distribution in different exposure categories increase the power of BF, SVLM, and DRM. Our results highlight vQTL detection methods that perform optimally under realistic simulation settings and show that their relative performance depends on the phenotype distribution, allele frequency, sample size, and the type of exposure in the interaction model underlying the vQTL.

Genome dosage alteration caused by chromosome pyramiding and shuffling effects on karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids.

KEY MESSAGE: We developed an array of Zea-Tripsacum tri-hybrid allopolyploids with multiple ploidies. We unveiled that changes in genome dosage due to the chromosomes pyramiding and shuffling of three species effects karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids. Polyploidy, or whole genome duplication, has played a major role in evolution and speciation. The genomic consequences of polyploidy have been extensively studied in many plants; however, the extent of chromosomal variation, genome dosage, phenotypic diversity, and heterosis in allopolyploids derived from multiple species remains largely unknown. To address this question, we synthesized an allohexaploid involving Zea mays, Tripsacum dactyloides, and Z. perennis by chromosomal pyramiding. Subsequently, an allooctoploid and an allopentaploid were obtained by hybridization of the allohexaploid with Z. perennis. Moreover, we constructed three populations with different ploidy by chromosomal shuffling (allopentaploid × Z. perennis, allohexaploid × Z. perennis, and allooctoploid × Z. perennis). We have observed 3 types of sexual reproductive modes and 2 types of asexual reproduction modes in the tri-species hybrids, including 2n gamete fusion (2n + n), haploid gamete fusion (n + n), polyspermy fertilization (n + n + n) or 2n gamete fusion (n + 2n), haploid gametophyte apomixis, and asexual reproduction. The tri-hybrids library presents extremely rich karyotype heterogeneity. Chromosomal compensation appears to exist between maize and Z. perennis. A rise in the ploidy of the trihybrids was linked to a higher frequency of chromosomal translocation. Variation in the degree of phenotypic diversity observed in different segregating populations suggested that genome dosage effects phenotypic manifestation. These findings not only broaden our understanding of the mechanisms of polyploid formation and reproductive diversity but also provide a novel insight into genome pyramiding and shuffling driven genome dosage effects and phenotypic diversity.

Megaproject will chart human immune diversity.

Human Immunome Project's survey aims to improve drugs and vaccines.

Macular neovascularization and polypoidal choroidal vasculopathy: phenotypic variations, pathogenic mechanisms and implications in management.

Advances in imaging have led to improved ability to characterize variations in clinical sub-phenotypes of macular neovascularization (MNV) in Age-related macular degeneration (AMD). Polypoidal choroidal vasculopathy (PCV) was initially described based on characteristic features observed in indocyanine green angiography (ICGA) and was thought to be a distinct entity from AMD. However, subsequent careful observations based on confocal scanning laser ophthalmoscopy-based ICGA, optical coherence tomography (OCT) and OCT angiography have led researchers to appreciate similarities between PCV lesion and type 1 MNV in typical neovascular AMD. Concurrently, clinical trials have shown that anti-VEGF monotherapy can achieve favourable visual outcome in the majority of eyes with PCV. These learnings have led to a shift in the way PCV is managed over the past decade. Recent studies have supported the use of non-ICGA based imaging modality to screen for PCV and the adoption of anti-VEGF monotherapy as initial therapy for PCV. A focus of recent research has been in the understanding of the role of choroidal alterations in the pathogenesis of PCV. The concept of pachychoroid in leading to outer retinal ischemia has garnered increasing support. Future research in this area should evaluate the potential of choroidal morphology in guiding personalized therapy in PCV.

Variable Phenotype of Congenital Corneal Opacities in Biallelic CYP1B1 Pathogenic Variants.

PURPOSE: The aim of this study is to describe the variable phenotype of congenital corneal opacities occurring in patients with biallelic CYP1B1 pathogenic variants. METHODS: A retrospective chart review was conducted to identify patients with congenital corneal opacities and CYP1B1 pathogenic variants seen at UPMC Children's Hospital of Pittsburgh. Ophthalmic examination, high-frequency ultrasound, anterior segment optical coherence tomography, histopathologic images, and details of genetic testing were reviewed. RESULTS: Three children were identified. All presented with raised intraocular pressure. Two patients showed bilateral limbus-to-limbus avascular corneal opacification that did not resolve with intraocular pressure control; 1 showed unilateral avascular corneal opacity with a crescent of clear cornea, iridocorneal adhesions, iridolenticular adhesions, and classical features of congenital glaucoma in the fellow eye (enlarged corneal diameter, Haab striae, and clearing of the corneal clouding with appropriate intraocular pressure control). The first 2 patients were visually rehabilitated with penetrating keratoplasty. Histopathology revealed distinct features: a variably keratinized epithelium; a thick but discontinuous Bowman-like layer with areas of disruption and abnormal cellularity; Descemet membrane, when observed, showed reduced endothelial cells; and no pathological changes of Haab striae were identified. Two patients had compound heterozygous pathogenic variants in CYP1B1 causing premature stop codons, whereas 1 was homozygous for a pathogenic missense variant. CONCLUSIONS: Congenital corneal opacities seen in biallelic CYP1B1 pathogenic variants have a variable phenotype. One is that commonly termed as Peters anomaly type 1 (with iridocorneal adhesions, with or without iridolenticular adhesions) and the other is a limbus-to-limbus opacity, termed CYP1B1 cytopathy. Clinicians should be aware of this phenotypic variability.

Stratification of the risk of ovarian dysfunction by studying the complexity of intermediate and premutation alleles of the FMR1 gene.

FMR1 premutation female carriers are at risk of developing premature/primary ovarian insufficiency (POI) with an incomplete penetrance. In this study, we determined the CGG repeat size among 1095 women with diminished ovarian reserve (DOR) / POI and characterized the CGG/AGG substructure in 44 women carrying an abnormal FMR1 repeat expansion number, compared to a group of 25 pregnant women carrying an abnormal FMR1 CGG repeat size. Allelic complexity scores of the FMR1 gene were calculated and compared between the two groups. In the DOR/POI cohort, 2.1% of women presented with an intermediate repeat size and 1.9% with a premutation. Our results suggest that the risk of POI is highest in the mid-range of CGG repeats. We observed that the allelic score is significantly higher in POI women compared to the pregnant women group (p-value = 0.02). We suggest that a high allelic score due to more than 2 AGG interspersions in the context of an intermediate number of repetitions could favor POI. Larger studies are still needed to evaluate the relevance of this new tool for the determination of the individual risk of developing POI in women with abnormal number of CGG repeats.

Digenic inheritance accounts for phenotypic variability in amelogenesis imperfecta.

Amelogenesis imperfecta (AI) represents a group of clinically and genetically heterogeneous disorders that affect enamel formation and mineralization. Although AI is commonly considered a monogenic disorder, digenic inheritance is rarely reported. In this study, we recruited two nonconsanguineous Chinese families exhibiting diverse phenotypes of enamel defects among affected family members. Digenic variants were discovered in both probands. In family 1, the proband inherited a paternal frameshift variant in LAMA3 (NM_198129.4:c.3712dup) and a maternal deletion encompassing the entire AMELX gene. This resulted in a combined hypoplastic and hypomineralized AI phenotype, which was distinct from the parents' manifestations. In family 2, whole-exome sequencing analysis revealed the proband carried a maternal heterozygous splicing variant in COL17A1 (NC_000010.11 (NM_000494.3): c.4156 + 2dup) and compound heterozygous variants in RELT (paternal: NM_032871.4:c.260A > T; maternal: NM_032871.4:c.521 T > G). These genetic changes caused the abundant irregular enamel defects observed in the proband, whereas other affected family members carrying heterozygous variants in both COL17A1 and RELT displayed only horizontal grooves as their phenotype. The pathogenicity of the novel COL17A1 splice site variant was confirmed through RT-PCR and minigene assay. This study enhances our understanding by highlighting the potential association between the co-occurrence of variants in two genes and variable phenotypes observed in AI patients.

Phenotypic variability in Joubert syndrome is partially explained by ciliary pathophysiology.

INTRODUCTION: Joubert syndrome (JS) arises from defects of primary cilia resulting in potential malformations of the brain, kidneys, eyes, liver, and limbs. Several of the 35+ genes associated with JS have recognized genotype/phenotype correlations, but most genes have not had enough reported individuals to draw meaningful conclusions. METHODS: A PubMed literature review identified 688 individuals with JS across 32 genes and 112 publications to bolster known genotype/phenotype relationships and identify new correlations. All included patients had the "molar tooth sign" and a confirmed genetic diagnosis. Individuals were categorized by age, ethnicity, sex and the presence of developmental disability/intellectual disability, hypotonia, abnormal eye movements, ataxia, visual impairment, renal impairment, polydactyly, and liver abnormalities. RESULTS: Most genes demonstrated unique phenotypic profiles. Grouping proteins based on physiologic interactions established stronger phenotypic relationships that reflect known ciliary pathophysiology. Age-stratified data demonstrated that end-organ disease is progressive in JS. Most genes demonstrated a significant skew towards having variants with either residual protein function or no residual protein function. CONCLUSION: This cohort demonstrates that clinically meaningful genotype/phenotype relationships exist within most JS-related genes and can be referenced to allow for more personalized clinical care.